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I .  INTRODUCTION 


This  report  describes  a  model  for  assessing  the  effects  of  scat¬ 
tering  and  absorption  of  infrared  radiation  by  atmospheric  aerosols  and 
gases  on  the  ranges  at  which  surface  targets  can  be  detected  against 


their  natural  backgrounds  by  a  surface-based  FLIR  system.  In  partic¬ 
ular,  equations  are  fitted  to  curves  of  the  8-12  um  band  averaged 


molecular  (range  dependent)  absorption  and  aerosol  extinction  coeffi- 

1  2 

cients  calculated  using  LOWTRAN  5  and  the  Navy  aerosol  model  (to 


appear  in  LOWTRAN  6),  respectively,  for  differing  meteorological  condi¬ 
tions.  These  equations  are  implemented  via  a  Hewlett-Packard  model  9845 
computer  program  ( listed  in  Appendix  form)  which  calculates  the  detec¬ 
tion  ranges  of  targets,  using  as  inputs  the  FLIR  system  parameters  and 
measurable  meteorological  parameters. 


II.  DETECTION  MODEL 


The  FLIR  system  minimum  detectable  temperature  difference  at  a 
spatial  frequency  v,  MDT(v),  is  that  temperature  difference  between  a 
square  target  of  angular  subtense  (2v)  ^  and  a  uniform  background  which 


1.  F.X.  Kneizys,  E.P.  Shettle ,  W.O.  Gallery,  J.H.  Chetwynd,  Jr.,  L.W. 
Abreu,  J.E.A.  Selby,  R.W.  Fenn  and  R.A.  McClatchey,  "Atmospheric 
Transmittance/Radiance :  Computer  Code  LOWTRAN  5,"  AFGL-TR-80-0067, 

21  February  1980. 

2.  Private  communications,  Stuart  G.  Gathman,  Code  4320,  Naval  Research 
Laboratory. 
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corresponds  to  a  50%  probability  of  detection  for  the  average  observer* 
It  defines  the  system's  ability  to  discern  an  object  which  may  be  too 
small  to  spatially  resolve  but  is  hot  enough  to  exceed  the  detection 
threshold  for  one  detector. 

MDT(v)  can  be  calculated  from  am  expression  given  by  Mash  and  Wood 
as 


MDT(v)  =  4v ( SNRT) ( NET ) 


(1) 


where 


SNRT  =  perceived  signal-to-noise  threshold  for  detection, 
(usually  set  to  2.5) 

NET  =  noise  equivalent  temperature  CC) 

r  =  system  resolution  (mrad) 

s 

Ax  =  detector  in-scan  angular  subtense  (mrad) 

Ay  =  detector  cross-scan  angular  subtense  (mrad) 
t@  =  eye  integration  time  (usually  set  to  0.1  sec) 

F^  =  Frame  rate  (frames/sec) 


Equation  ( 1 )  is  applicable  to  common-module  FLIR  systems  and  when  the 
system's  modulation  transfer  function  is  due  only  to  the  detector  and 
optics. 


3.  C.  Nash  and  J.T.  Wood,  "Calculation  of  MRT  and  MDT  for  Thermal 

Imaging  Sensors,"  DCS  Corporation  Report  to  NSWC  under  Contract  No. 
N6092 1-81 -M-49 1 7 ,  24  August  1981. 
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Real  targets  are  generally  complex  In  shape.  Here,  for  modeling 


purposes,  the  target's  composite  areas  are  sunned  to  give  an  equivalent 

rectangular  area  of  width  and  height  Hfc.  The  target  solid  angle 

2 

subtense  is  found  by  dividing  its  area  by  its  range  squared  (R  ),  and 
the  angular  subtense  8  is  the  square  root  of  this.  The  target  spatial 
frequency  is  then  given  by  v  =  (20)  \  Equation  (1)  can  now  be  written 


MDT(v)  =  C1  v 


(,  ♦  c,.*)'" 


MDT(R)  =  C 


1  2<HtWt)V2 


1  +  C 


2  4HtWt 


1/2 


(2a) 


(2b) 


where 


and  C2 


are  system  constants  given  by 


17. 84 (NET) 


(3a) 


(3b) 


The  perceived  temperature  difference,  AT^ ( R ) ,  between  a  target  at 
range  R  and  its  background  is  the  actual  temperature  difference  At 

A 

degraded  by  the  atmospheric  transmittance,  <x(R)>,  band-averaged  over 
the  response  of  the  the  system,  ie. 


AT  (R)  -  AT  <T(R)> 
P  a 


(4) 


i 
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The  band-averaged  transmittance  is  defined  as 


<t(R)> 


-<B(R)>R 

e 


(5) 


where 


<B>  =  <$  >  +  <A  ( R)>  (6) 

a  m 


and  <8  >  and  <A  (R)>  are  the  band-averaged  aerosol  extinction  and  range- 
a  m 

dependent  molecular  absorption  coefficients  which  will  be  derived  in  the 
following  sections. 


The  detection  range  is  then  defined  as  that  range  where  AT  (R)  5 

P 

MDT(R)  or 


-<B(R)>R  „ 

AT  e  =  C, 

a  1 


2(HtV 


1/2 


1  + 


C2R 

4Ht\ 


1/2 


(7) 


Given  the  atmospheric  parameters,  target  dimensions  and  actual  tempera¬ 
ture  difference,  the  detection  range  is  determined  from  Eq  (7)  using  a 
Hewlett-Packard  model  9845  computer.  Basically,  the  computer  program 
(listed  in  the  Appendix)  iterates  the  difference  between  the  left  and 

right  sides  of  Eq  (7)  by  varying  R  until  the  difference  is  postive  and 
-5 


less  than  10 


In  the  iterative  process,  the  perceived  range  is  corrected  for 

blocking  by  the  horizon.  If  the  incremented  range  is  less  than  the 

optical  horizon,  R  ■  3.85  /H~ ,  where  H  is  the  sensor  height  in  metres, 
o  s  s 

there  is  no  blocking.  If  R  >  Rq  the  target  height  used  in  the  iteration 
process  is 


(8) 


III.  AEROSOL  EXTINCTION  MODEL 

The  aerosol  extinction  model  incorporated  into  SUBFLR  was  derived 

2 

using  the  Navy  aerosol  size  distribution  model  developed  at  the  Naval 
Research  Laboratory  for  inclusion  into  LOWTRAN  6.  'Hie  size  distribution 
model  (developed  from  empirical  fits  to  a  large  data  set  of  size  dis¬ 
tribution  measurements  in  a  variety  of  marine  environments)  is  the  sum 
of  three  log-normal  distributions.  Component  one  represents  continental 


aerosols,  and  components  two  and  three  represent  equilibrium  sea  spray 


particles  generated  by  the  24-hour  and  instantaneous  surface  wind 
speeds*  respectively,  given  by: 


n(r) 


(10) 


where  A^  *  2000  (AM)^ 

A_  =  5.866  (W  -  2.2) 

2  o 

A3  =  .01527  (W  -  2.2) 


(AM)  is  an  air  mass  parameter  varying  from  a  value  of  1.0  for  open  ocean 
to  a  value  of  10  in  coastal  areas.  Wq  and  W'  are  the  averaged  24-hour 
and  instantaneous  wind  speeds,  respectively.  If  Wq  and  W  are  unknown 
the  following  zonal  default  values  apply  with  Wq  =  W : 


SL  • 

Tropical 

=  4.1 

m/s 

b. 

Midlatitude 

summer 

=  4. 1 

m/s 

c  • 

Midlatitude 

winter 

=  10.3 

m/s 

d. 

Subarctic 

summer 

=  6.7 

m/s 

e. 

Subarctic 

winter 

=  12.4 

m/s 

Other  default  values  are:  if  5.866(W  -2.2)  <  0.5  then  A_  =  0.5,  and  if 

o  i 

0.01527(W'-2.2)  <  1.5  x  10-5,  then  A3  =  1.5  x  10-5. 


6 


In  Bq  (10)  ,  the  modal  radius  for  each  distribution  component,  is 


allowed  to  grow  with  relative  humidity  (RH)  according  to 


f  ,  f  <2-7  M/L00)-,  1V3  (11 

l  6(1  -  RH/ 100)  J  V 

The  contribution  to  the  total  extinction  by  each  component  can  be  writ¬ 
ten  as: 


B  (X). 
a  l 


f 0  ( X ,r ,m)  exp 
/ 

L  V  x/J 

r2  dr 


(12) 


where 


C, 

i 


0.001*  . 
i  Ai 


The  factor  f  1  in  the  expression  for  Ci  insures  a  constant  total  number 

of  particles  as  the  humidity  increases*  Qe(X,r,e)  is  the  cross  section 

for  extinction  normalized  to  the  spherical  particle  geometrical  cross 

section,  and  m  is  the  complex  refractive  index  at  a  wavelength  X.  The 

4  5 

refractive  index  is  allowed  to  change  from  that  of  a  dry  sea  salt  as 
the  particle  deliquesces  with  increasing  humidity. 

4.  G.  Hanel,  "New  Results  Concerning  the  Dependence  of  Visibility  on 
Relative  Humidity  and  Their  Significance  in  a  Model  for  Visibility 
Forecast,"  Beitrage  zur  Physik  der  Atmosphere,  44,  137-167,  1971. 


5.  F.E.  Volz,  "Infrared  Refractive  Index  of  Atmospheric  Aerosol 
Substances,"  Appl.  Opt.,  _1_1,  755-759,  1972. 


The  Navy  model  provides  precalculated  values  of  the  parameters 
6a(X)i/Ci  at  discrete  wavelengths  for  relative  humidities  of  50,  85,  95 
and  99%  from  which  the  average  extinction  coefficient 


<B  >. 

a  i 


h  /■ 


B  U).  dX 
a  i 


(13) 


for  a  wavelength  band  AX  can  be  determined. 

Shown  in  Fig  1  are  the  values  of  <B  >./C.  calculated  from  the 
’  ail 

discrete  wavelength  values  provided  by  the  Navy  model  at  the  four  rela¬ 
tive  humidities.  Also  shown  are  the  least  squares  fit  to  the  data 
points  by  the  equation 


<Ba>i/Ci  100  -  RH 


+  b. 

i 


(  14) 


-ne  values  of  a^  and  b^  determined  from  the  regression  analyses  allow 
the  8-12  ym  aerosol  extinction  coefficient  to  be  written  as: 


t 


x  10-9) 


(S 


-3 

5.86  (W  -  2.2)  |  —  +  3.69  x  10  '  |  + 

o  l  100  -  RH 
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RELATIVE  HUMIDITY  (%) 


Figure  1.  8-12  u m  band  averaged  values  of  <Ba>VC 

calculated  from  the  Navy  model  at  relative  humidi¬ 
ties  of  50,  80,  95  and  99%.  Connecting  lines 
indicated  least  squares  fit  to  the  data  point. 


0.01527  (W  -  2.2) 


(iS“= 


km 


-1 


(15) 


In  a  similar  manner  a  relationship  for  the  extinction  coefficient  at 
0.55  |<n  is  determined  (see  computer  listing  line  1040).  If  the  visual 
range  (visibility)  Ry  is  known,  the  aerosol  extinction  coefficient  can 
be  normalized  by  the  scaling  factors 


V 


sf 


3.912 

R  0  (0.55  pm) 
v  a 


(16) 


This  normalization  essentially  adjusts  the  total  number  of  particles  to 
match  the  observed  visibility.  The  above  relationships  apply  to  rela¬ 
tive  humidities  between  50  and  99%  and  assume  horizontal  homogeneity. 

If  the  humidity  is  unknown  a  default  value  of  80%  should  be  used. 


IV.  RANGE-DEPENDENT  MOLECULAR  ABSORPTION  COEFFICIENT 


In  addition  to  water  vapor  absorption  (band  and  continuum),  molecu¬ 
lar  absorption  in  the  8-12  pm  band  is  influenced  by  uniformly  mixed 
atmospheric  gases  (C02,  N2,  CH^,  CO,  02»  and  O^).  If  the  absorber 
quantities  along  the  path  are  known,  the  atmospheric  transmittance  can 
be  determined  for  the  path  from  the  spectral  absorption  properties  of 
the  various  absorbers.  The  total  broadband  transmittance  is  calculated 


10 


by  multiplying  the  individual  transmittances  of  each  contributing  ab¬ 
sorber  averaged  over  a  narrow  spectral  band,  61,  and  integrating  these 
composite  results  over  the  broad  spectral  band  AX*  Thus 


<t  (R)>  =7-  A*.  T.(R,«M]dX  (17) 

®  AX  J  3  3 

X 

The  molecular  absorption  model  employed  in  SUBFLR  has  been  derived  from 
the  LOWTRAN  5  program  of  Kneizys  et  al1.  The  8-12  um  band  averaged 
absorption  coefficients  computed  from  <xm(R)>  of  Eq  (17)  according  to 
the  relationship 


<A  ( R)>  =  -  7  In [<T  ( R) > ]  (18) 

m  R  in 

are  shown  in  Fig  2  as  functions  of  range  for  differing  absolute  humidi¬ 
ties.  The  curves  of  constant  absolute  humidity  were  computed  for  a 
surface  level  (constant  pressure)  path  using  several  different  tempera¬ 
tures  and  relative  humidities.  The  contributions  to  transmittance  by 
the  uniformly  mixed  gases  were  accounted  for  by  assuming  absorber 
amounts  equivalent  to  the  1962  Standard  Atmosphere.  The  circled  curves 
indicate  different  absorptions  calculated  with  the  same  absolute  hu¬ 
midity  which  resulted  from  different  combinations  of  temperature  and 
relative  humidity.  These  differences  in  absorption  are  most  likely 
caused  by  how  the  band  and  continuum  contributions  depend  upon  the 
actual  temperature  and  relative  humidity. 


i 
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Molecular  Absorption  Coefficient  (km 


Figure  2.  8-12  yin  band  averaged  absorption  coefficient 

versus  range  for  different  absolute  humidities. 
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In  attempting  to  obtain  a  functional  relationship  for  the  absorp¬ 
tion  curves  in  Fig  2,  it  was  found  that  the  residuals  of  the  least 
squares  fit  to  the  data  were  a  minimum  with  the  form 

<A  (R)>  -  a  +  a,  log  R  +  a  ( log  R)2  (19) 

m  o  1  2 

where  the  coefficients  oqI  and  43:6  functions  of  absolute  humidity. 
Functional  forms  for  the  band-averaged  absorption  coefficient  at  speci¬ 
fic  ranges  can  also  be  derived,  frcm  which  values  of  a  ,  t.,  and  a.  can 

o  i  3 

be  determined  for  different  absolute  humidities.  In  Fig  3,  the  band- 
averaged  absorption  coefficients  calculated  at  1,  10  and  100  km  are 
shown  as  a  function  of  absolute  humidity.  The  absorption  coefficient  at 
each  range  can  now  be  represented  in  terms  of  absolute  humidity.  Ah,  by 
a  second-order  polynomial  equation  of  the  form 

B(R)  -  a  +  a, (Ah)  +  a, (Ah) 2  (20) 

o  1  2 

where  the  coefficients  determined  for  the  three  ranges  are  given  in 
Table  1 . 


Table  1.  Coefficients  for  three  ranges. 


1  8.0024E-2  1 . 2360E-2  5.3855E-4 

10  1 . 5916E-2  7 . 8820E-2  4.9786E-4 

100  4. 8333E-3  4.7687E-3  4.8820E-4 
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Molecular  Absorption  Coefficient  (km 


Figure  3.  8-12  un  band  averaged  molecular  absorption 

coefficient  versus  absolute  humidity  at  ranges  of  1, 
10  and  100  km. 
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Equating  Eq  (19)  and  (20)  at  each  of  the  three  ranges,  the  coefficients 
in  Eq  (19)  can  be  determined  as 


®c  *  B1  (20a) 

*1  “  B2  ~  ao  "  a2  (20b) 

a2  =  B3  -  2B2  +  B1  (20c) 

where  B^ ,  B„,  B^  correspond  to  the  values  of  B(R)  at  1,  10  and  100  km, 
respectively. 


V.  SAMPLE  OUTPUT 


After  the  program  is  loaded  into  memory,  the  HP-9845  computer 
interactively  queries  the  user  for  the  required  inputs.  As  the  asked- 
for  inputs  are  typed  in,  the  "CONTINUE"  key  is  pressed  for  sequencing  to 
the  next  input  and  initiating  the  range  calculation  and  display/print¬ 
out.  A  sample  printout  for  a  current-generation  PLIR  operating  against 
a  large  cruiser  in  the  open  ocean  is  shown  below.  The  meteorological 
parameters,  temperature  difference  and  sensor  height  were  arbitrarily 

chosen  and  resulted  in  a  50%  probability  detection  range  of  14.93  km. 

In  this  case  the  horizon  obscured  4.6  m  of  the  target's  19-m  height. 


(sample  computer  printout) 


NEAR  SURFACE,  HORIZONTAL  DETECTION  RANGE  CALCULATION 


HM  =  1 

Nirid  24  Hr  A"e.=  4.1  MS  Wind  Inst.  =  4.1  MS  V  l  s  i  b  l  1  1 1  ;.■  =  1 5 . 0  km. 

Delta  T  =  3.0  C  FLIR  Ht=  3.0  Metres  CONST.  1  =2.036E-01  CONST.  2  »1 . 568E-01 
P h  =  80. 0  Temp. =16.0  Deg.  C  Ah=10.9  Gm  M  3  Aero  E<t=  .170 
Target  Heights  19.00  Width=160.00  Metres 

TOTAL  EXT=4 . 60 1  RANGE=  14.4  KM  OP  7.8  NM 

Obscurations  4.1  Metres  Effective  Target  Height =14. 93  Metres 
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10  OPTION  BASE  1  ‘Name**SUBFLR***15  SEPT.  1982  **MRP*» 

20  !  MODIFIED  TO  CfiLC  EXT.  IN  SUBROUTINE  "Hughes"  31  Rug.  1982 

30  ,  - 

40  !  **This  program  requires  the  input  of  the  various  FLIR  parameters  as  uell 

50  !  as  the  t arget -t o-background  temperature  difference.  Designed  for  surface 

60  !  detection  ranges.  Allows  for  target  obscuration  for  beyond  the  horizon 

70  !  condi t i ons. 

80  i  - 

90  PRINTER  IS  0 

100  PRINT  CHRS ( 27 ' & “ &  1 20T " 

110  PRINTER  IS  16 

120  INPUT  "PRINTER  SELECTION.  P  for  printer,  S  for  screen", X» 

130  IF  X*»"P"  THEN  PRINTER  IS  0 
140  IF  X*=“S"  THEN  PRINTER  IS  16 

150  PRINT  “  - 


160  PRINT  TAB* l 1>; "NEAR  SURFACE,  HOR I 20NTRL  DETECTION  RANGE  CALCULATION" 
170  PRINT  "  - 


180  Start:  i  ♦+♦♦♦♦♦+♦+♦♦♦♦♦♦♦♦ 

190  RAD 

200  i  ♦♦♦+♦♦♦♦♦♦+♦♦+♦♦+♦♦♦ 

210  Hrange:  1  *  Calculates  horizontal  detection  range  from  FLIR  to  target.  Uses 
range  dependent  molecular  extinction. 

220  1  ******** 

230  INPUT  "DELTA  T>C>,  SENSOR  HE  I GHT < Met  res > “ , Td i ff, Hsensor 

240  INPUT  "SYSTEM  CONSTANTS  C 1 , C2 " , Cs 1 , Cs2 

250  INPUT  "RH, TEMP. < DEG.  C<",Px, Temp 

260  GOSUB  Hughes 

270  Aerl=Betair 

280  PRINT  USING  290; Td i ff , Hsensor , Cs 1 , Cs2 

290  IMAGE  2X, "Delta  T  =  ".DD.D,"  C  FLIP  Ht  =  " , DDD. D, "  Metres  CONST.  1  =" , D. DDDE, 

"  CONST.  2  , D. DDDE 

380  GOSUB  Abshum 

310  PRINT  USING  320; Px , Temp, Ah, Aer 1 

320  IMAGE  2X, "Ph  =  " , DD. D,  “  Temp. « ” , DD . D, "  Deg.  C  Ah«",DD.D,"  G» 'M  3  Aero  Ext»",D 
.DDD,"  Mol.  Ext . *"D. DDD 

330  INPUT  "Target  Height, Width  in  me* res " , Ht ar get , U 
340  PRINT  USING  350; Ht arget  ,  W 

350  IMAGE  2X, "Target  He i ght « " , DDD. DD, "  Wi dt h=" , DDD. DD"  Metres" 

360  R*0 

370  Rd»l 

380  i  - 

390  Loop:  R»R+Pd 

400  i  - - - - 

410  Hobsc»< P-3. 85*S0PtHsensor )  ■  2  14.8225 

420  IF  P< 3. 85*SQR' Hsensor >  THEN  Hobsc»0 

430  H»Ht arget -Hobsc 

440  DISP  . Hobsc, H 

450  GOSUB  Rmolec 

460  T«Aerl+Amol 

470  IF  H* «0  THEN  0"erhoriz 

480  Al«Tdi f f *EXPc -T*P> 

490  Hew”Rr SQR<  H*M )'2 

500  A2«Csl*New*S0Ptl*Cs2*Neu'  2  • 

510  A-A1-A2 

520  IF  A*0  THEN  GOSUB  Pdel 

530  IF  <A>0>  AND  < A'. 00001*  THEN  GOTO  Prntout 
540  GOTO  Loop 

550  PRINT  TAB<  3 ) , . - - - - - - - 


560  GOTO  Start 
570  '  *»**•••******» 

580  Rde 1 :  '  • 

590  ■  ************** 

600  R-R-Rd 
610  Rd>Rd'10 
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628  RETURN 
630  !  ********* 

640  Abshum  I  1  * 

650  I  ********* 

660  0*18x8.315 

670  Esat  *6.  105*EXP(  25. 22*  ( Temp  i  Terup+273.  2  '  >-5. 31  *L0G'  1  +  Temp-' 273.  2  >  ) 

680  Abhum*Rx*E*at  *0  - >  273 . 2  +  Temp  > 

688  Ah*Abhum 

700  RETURN 

710  !  ********** 

720  Prnt  out :  •  * 

730  1  ********** 

740  Xt  t 1 *R*T 
750  Rnm  =  R  1 . 853 
760  PRINT 

770  PRINT  USING  780 ;  Xt  t  1  ,  R ,  Rr.m 

780  IMAGE  2X, “TOTAL  EXT* “ , D . DDD . "  RANGE* ", DDB. D. "  HI  OP  " , HDD .  Il , "  NM" 

730  PRINT  USING  800 1  Hobs c , H 

800  IMAGE  2X,  “  Obsc  ur  at  t  on*  "  .  BDD.  B,  “  Metres  Effect  lue  Target  He  t  ght  =  "  ,  BB .  BB,  "  Me 
t  re  »  " 

810  GOTO  550 
830  Rmotec:  1 

850  B 1 =8 . O024E-2* 1 . 2 36E -2* Ah  +  5 . 3855E-4*Ah* Ah 

360  B2  = 1 . 531 6E-2*7. S32E- 3* Ah* 4. 3786E -  4* Ah*Ah 

870  B  3  =  4 . 833  3E-  3+4 . 7687E-3*Ah+4 . 882C'E-4*Ah*Ah 

880  Alpha0=61 

830  A 1 pha2  = < B3-2+B2+B l >.  2 

300  A 1 phal =B2-A 1 ph a0-A ! pha2 

310  Amo  1 =A I ph a0  +  Al phal *LGT • P j+A 1 pha2*LGT ' P •  2 
320  RETURN 

330  Hughes:  1  Calculates  aerosol  extinction*  per  NRL  equations.  31  Aug.  1382 
340  INPUT  "AM  <1-10  •".An. 

350  PRINT  "  AM=“;Am 

360  R<l=100-P 

370  INPUT  "24  Hour  Wind  speed  a"e.  Metres  sec.",W0 
380  INPUT  "Current  measured  wind  speed  Metres  ssc.",Ws 
330  02*5.866*040-2.2.) 

1000  C3* . 0 1 527*t Ws-2. 2  > 

1010  IF  C2  .5  THEN  C 2 * .  5 

1020  IF  C30.4E-5  THEN  C3-1.4E-5 

1030  i  - CALC.  EXTINCTION  AT  0.55  - 

1040  Be»  a55= • 4 . 74E-2  R  1 +2. 3E-4  ■  * Am * Am +  C2*'  1 . 48E-2  R  a 1  +  3. 08E-4  > +C  3*' 7. 73  R • 1 ♦ . 42 
8  > 

1050  INPUT  "Visual  Range  in  In.  it'  unk  noun  input  -l“,Rv 
1060  IF  P"  0  THEN  Vsf *3. 312  •Pv*Beta55> 

1070  IF  Ru  0  THEN  Vsf  *  1 

1080  i  --CALC.  IR  EXTINCTION  - 

1030  Al*'  1.1 38E-3  R  1  +  3. 32E-6  '»Am*Am 
1100  A2=  <"  4 .  1  3E  -  3  P-  1  +  3. 63E-5  >  *C  2 
1110  A3*<  8.  308  RO  +  .431>*C3 
1 120  Bet  ai  r  *•:  Al  +  A2  +  A3  )*V*f 
1130  PRINT  USING  1 1 40 ; W0 , Ws , Pu 

1140  IMAGE  2X,“Wind  24  Hr  Ave . = " , DP . B, "  MS  Wind  I ni t . * " , DB . D, "  M  S  V l s i b i 1 l t u  = 
" , DB . D, "  km." 

1150  RETURN 

1160  Overhortz:'  *********** 

1170  PRINT 

1180  PRINT  "TARGET  IS  OVER  THE  HORIZON,  OUT  OF  VIEW" 

1130  PRINT  "  - - 


1200  GOTO  Start 
1210  END 
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